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Charmonium states can decay into pairs of D and D mesons if their masses are above the al- 
lowed decay thresholds. In general cc states near threshold will also undergo mixing with DD 
molecular (or tetraquark) states, by creation and annihilation of light quark-antiquark pairs. The 
investigation of such effects sheds light on the higher Fock state contributions to charmonium 
wavefunctions and on mass shifts, relative to a scenario where such mixing effects are neglected. 
A variational approach is applied to a mixing matrix between operators of both sectors, of cc and 
of DD molecular type. The efficient calculation of several diagrams appearing in this matrix re- 
quires all-to-all propagators, which are realized by sophisticated stochastic estimator techniques. 
The runs are performed on np — 2 24 3 x 48 lattice volumes with m n m 280 MeV, using the non- 
perturbatively improved clover Wilson action, both for valence and for sea quarks. 
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1. Introduction 

Nonperturbative simulations have demonstrated that quantum chromodynamics (QCD) pre- 
dicts the confinement of objects that carry a color charge. The hadronic color singlet states with 
baryon number B = are called mesons. Within the nonrelativistic quark model these contain 
exactly one quark and one antiquark. However, at least in principle, QCD offers the possibility to 
construct such states entirely out of gluons (glueballs), out of quark, antiquark and gluons (hybrids), 
out of two quarks and two antiquarks (molecules or tetraquarks) and even higher Fock states. In 
general these states will undergo mixing and not all these states will be stable against strong decays. 

The charmonium sector provides us with a particularly rich laboratory for the study of these 
effects. Several experimentally discovered charmonium states are suspected to contain large con- 
tributions from non-cc configurations m, ^, |p. Here we wish to address the contribution of cqqc 
molecules 1 to S-wave charmonia and of cc states to DD bound states. 

QCD predicts the spectrum within a sector of fixed J PC , isospin /, strangeness S and charm C. 
Classifying the resulting states according to their (presumed) quark content will be model depen- 
dent. However, there are examples where clearly one Fock state dominates over another and where 
mixing is small. For instance the (J PC = 0~ + , / = 5 = C = 0) r] c state is clearly different from the 
r\' state with exactly the same quantum numbers and in this case dynamic mixing effects turn out 
to be negligible such that we can identify the leading Fock contribution to the former state as 
cc while the 7]' is composed of strange and light quark-antiquark pairs. 

With such caveats in mind, we will somewhat carelessly use the term "T7 c -DiD*-mixing" 
when we mean mixing between states that couple to c/sc-type operators and states created by, 
{ c Y5Yi<i)(cYiq)- These "unperturbed" states can be accessed on the lattice by appropriate interpolat- 
ing fields. Their mixing is studied by applying a variational approach. One of the problems is that, 
unless the two sectors completely decouple from each other, in the limit of large Euclidean times 
the states created by all these operators will decay into the same ground state: the ground state 
created by the cc type operators, that we call the unmixed r\ c , will intrinsically already contain a 
D[D* contribution. However, if the perturbative approach that we outline below is justified, then 
this implicit mixing [Q] will first appear at second order in an expansion parameter A, while the 
transition matrix element between the two sectors (explicit mixing) will be of order X. Conse- 
quently, at intermediate Euclidean times implicit mixing might be negligible, reducing the model 
dependence of our ansatz. 

Our perturbative set-up is as follows. We expect the physical r] c wavefunction at first order in 
a parameter A to read, 



with a normalisation factor jV and a (hopefully) small 2 coupling constant A appearing in the 
mixing vertex of the Hamiltonian H\. While we do not know the functional form of H\ or of 
the unperturbed wavefunctions, we can evaluate all the relevant matrix elements on the lattice. 

'We use the term "molecule" synonymous also for tetraquarks. 

2 In some cases A might be large but we should be able to detect this through implicit mixing, i.e. decay of the state 
created by the cqqc operators into the cc states at small Euclidean times. 




(1.1) 
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One important thing to note here is the dependence of the mixing on the light quark mass. With 
decreasing m q , the denominator obviously becomes smaller, but the mixing matrix element in the 
numerator is expected to increase, since the probability for creating a light quark-antiquark pair 
should be inversely dependent on the light quark mass. Therefore we expect mixing effects to 
increase at smaller light quark masses. 




Figure 1: The mixing matrix. Solid lines represent charm quarks, wiggled lines light quarks. 



2. Simulation 

We aim to calculate the coefficients in the expansion of the QCD eigenstates into the trial 
interpolating fields nonperturbatively, by diagonalizing a matrix of cross correlators including both 
cc and molecular operators. For each type of operator we apply three different types of smearing: 
local, narrow and wide, where these terms indicate the number of Wuppertal smearing steps with 
8 = 0.3, employing spatial APE smeared (tiape = 15, a = 0.3) parallel transporters to smooth the 
trial wavef unctions. The number of fermion field smearing iterations is determined by optimizing 
the effective masses separately for each of the two sectors. 

In figure [I] we sketch the structure of the mixing matrix. The different smearing levels are 
omitted for the sake of clarity. Solid lines represent charm quark propagators and wiggly lines light 
quark propagators. The prefactors are due to the two mass degenerate light sea quark flavors. The 
upper left corner contains the cc, the lower right corner the molecular sector. Nonvanishing off- 
diagonal elements indicate mixing. The spatial separation within the molecular operators was tuned 
by maximizing the magnitude of the off-diagonal element. The optimal value was r = 4a « 0.3 fm. 

The charm-anticharm annihilation diagrams were omitted in the present study where we focus 
on cc-molecule mixing near threshold. We studied these previously in the context of r} c -T}' mixing 
where their effect turned out to be negligible [Q]. 

For the evaluation of the last two diagrams of the molecular sector, light all-to-all propagators 
are necessary. 0(~ 100) complex Z2 stochastic estimates per configuration were calculated for 
this purpose, with the application of sophisticated noise reduction methods like staggered-spin- 
color dilution and hopping parameter acceleration as already used in previous studies []6[ f|, ^]. 
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The generic form of our meson interpolators, centred around 
a position x, reads, 

M(x) = {cT M c) x , (2.1) 
while the molecular interpolators with separation r look like, 
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Y(x,r) = -L ((^c),(cr2 (7 ), +r + (-)'Xcr^),(^ c ) A . +r ) 



(2.2) 



Table 1: T structures of me- 
son and molecular interpolating 
fields. 



In table [l] we display the explicit T structures for the J PC = + , 1 and 1 ++ channels. 
The variational method consists of solving a generalized eigenvalue problem, 



C(t Q )- l / 2 C(t)C(t )- l / 2 V a =X a (t,t ) ¥ a , 



(2.3) 



see refs. 



10p for details. For sufficiently large times the eigenvalues and -vectors will approach 
their asymptotic values. The components of a given eigenvector can be interpreted as the coupling 
strengths of the corresponding operators to the state under consideration. 

Our strategy differs from that of many previous studies that utilized the variational method, 
in as far as our primary interest lies in the couplings and not only in the resulting spectrum. We 
first determine the eigenvalues of the three by three submatrices within each of the Fock sectors, 
separately, in order to obtain an "unperturbed" spectrum. This is then used to identify the affiliation 
of the eigenvectors to the eigenstates on the jacknifes samples. Finally, we will compare spectrum 
and eigenvector components with the mixing elements switched on to this unmixed reference point. 
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Figure 2: Mass spectra from the separate diagonalization of the submatrices within each sector. 

Our runs are performed on 100 effectively de-correlated np = 2 24 3 x 48 configurations gen- 
erated by the QCDSF collaboration [pd|], with a lattice spacing a 0.076 fm, obtained from the 
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chirally extrapolated nucleon mass. We employ the same non-perturbatively order a improved 
clover action for the valence quarks. The charm quark mass was set by tuning the spin-averaged 
charmonium mass \{m^ c +3m/m/) to its experimental value. The pion on these lattices is reason- 
ably light, m n « 280MeV, to embrace the above mentioned dependence of the mixing on the light 
quark mass. 

Computations were performed using the Chroma software library Jl2| ] partly on the local HPC 
cluster and partly on the BlueGene/P of the Julich Supercomputer Center. 



3. Spectrum 



An extra benefit bucked off by our analysis is the mass spectrum in the investigated channels. 
The separate diagonalization of the three by three submatrices provides us with at least four reliable 
eigenvalues, two for each subsector. However, since the molecular channels typically are rather 
noisy, we are only able to extract the ground states there, within reasonable errors. So we are left 
with three states in each channel, plotted in figure ||. For the molecular masses we give two data 
points: the left ones are from the diagonalisation procedure, the right ones represent the sums of 
the masses of the corresponding pairs of non-interacting D mesons. Note that in the 1 ++ channel, 
the radially excited % c \ is heavier than the molecular state, in contrast to the other channels. If we 
consider the fact that our pion is about 130-140 MeV too heavy, the mass of the molecular state is 
indeed consistent with the X (3872), which most likely has J PC 
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Figure 3: Effective masses of the eigenvalues of the submatrices (left hand side) and of the full matrix (right 
hand side) in the 1 channel. 



4. Mixing 

Equipped with the reference eigenvalues from the submatrices, we go for the diagonalization 
of the full six by six matrix. However, due to limited statistics, we find this to be numerically 
unstable and restrict ourselves to the sub-basis Mi oca i ,M narrow ,Yi oca i ,Y narrow . We discuss the vector 
state as one example. In the left of figure || we display the effective masses from the diagonalization 
of the two submatrices. The data points for J and its radial excitation are from the cc submatrix, 
the ones for the D{Dq from the molecular sector. 
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Figure 4: Eigenvector components from the diagonalization of the full matrix in the 1 channel. 



The unmixed reference points can also 
be found in black color in the right of fig- 
ure [| In addition, the two lowest lying ef- 
fective masses from the diagonalization of 
the full matrix are shown there. We are able 
to identify these two states with J and 
D\Dq, respectively. Interestingly, the J /W 
state is not found in the diagonalisation of 
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the four by four system. Table 2: Eigenvector components in the full basis. 

The corresponding eigenvector components are plotted in figure f| for the // \ff (left) and for the 
D]Do molecule. The J receives the dominant contribution from the local cc operator. However, 
the molecular configurations seem to contribute significantly too. The D\Dq state in contrast only 
contains small (but non-vanishing) cc admixtures. 

In table || we summarize the results for all channels that we investigated. In each of them we 
see significant mixing effects. The large molecular contribution to the Xcl is particularly noteworthy 
to mention. 

5. Conclusion & Outlook 

In all of the investigated channels (0~ + , 1 , and 1 ++ ) we detect significant mixing effects 
between cc and four-quark states. Although the precise values of the eigenvector components 
should not be taken too seriously since the operator basis is rather small and thus may miss non- 
negligible parts of the physical wavefunction, our analysis clearly substantiates the assumption of 
charmonium states having a rich Fock structure. 

In the near future this study will be extended to other channels including additional interpolat- 
ing fields, e.g. a molecular operator containing two mesons in a relative P-Wave. 
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